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INTRODUCTION 


Ionization  in  the  earth’s  atmosphere  is  produced  primarily  by  solar  illumination  and 
incoming  energetic  particles.  Ionization  produced  by  solar  illumination  is  well  understood 
and  can  be  accurately  modeled.  But  at  high  latitudes,  particularly  at  night,  the  dom¬ 
inant  source  of  ionization  may  be  precipitating  protons  and  electrons  with  energies  less 
than  about  100  keV.  This  makes  modeling  high  latitude  ionization  difficult  because  these 
particle  fluxes  may  be  temporally  and  spatially  varying.  Global  specification  of  the  distri¬ 
bution  of  ionization  at  high  latitudes  is  best  done  using  detectors  on  space-based  platforms 
capable  of  viewing  a  large  portion  of  the  earth  at  one  time.  This  method  involves  detect¬ 
ing  electromagnetic  radiation  originating  from  the  impact  of  charged  particles  into  the 
atmosphere.  This  radiation  may  be  in  the  X-ray,  ultraviolet,  or  visible  portion  of  the 
spectrum.  The  X-ray  and  ultraviolet  techniques  hold  great  promise  for  remotely  sensing 
ionization  because  of  the  absence  of  a  background  smission  source  that  prohibits  imaging 
of  aurorally  produced  emissions  in  the  sunlit  hemisphere.  However,  detailed  comparisons 
between  modeled  and  measured  auroral  emissions  are  required  to  assess  the  accuracy  of 
these  techniques. 

In  this  paper  we  compare  simultaneous  measurements  by  the  Chatanika  radar  and  the 
DMSP-F2  satellite  X-ray  detector  to  demonstrate  that  remote  X-ray  observations  can  be 
used  to  determine  the  ionization  produced  by  auroral  electron  precipitation.  The  accuracy 
of  the  X-ray  remote  sensing  is  assessed  by  a  quantitative  comparison  between  the  profiles  of 
ionization  measured  by  the  ground-based  incoherent-scatter  radar  and  the  profiles  derived 
from  the  satellite-based  X-ray  measurements.  Measurements  were  made  under  a  variety 
of  auroral  conditions,  including  dynamic  evening  aurora  and  sunlit  diffuse  aurora  in  the 
morning. 

The  next  section  of  this  paper  describes  the  method  by  which  auroral  electron  fluxes 
are  calculated  from  X-ray  measurements.  The  following  sections  describe  radar  and  satel¬ 
lite  data  obtained  during  two  auroral  events:  an  evening  pass  on  21  February  1979  and 
a  morning  pass  on  15  June  1978.  The  concluding  section  summarizes  the  results  and 
discusses  the  utility  of  the  X-ray  remote  sensing  method. 
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DETERMINATION  OF  AURORAL  ELECTRON  FLUXES  FROM 
X-RAYS 


Precipitating  electrons  interact  with  the  atmosphere  to  produce  ionization  and  various 
forms  of  electromagnetic  radiation.  Excitation  of  neutral  and  ionic  constituents  caused  by 
impact  of  primary  and  secondary  electrons  leads  to  optical  emissions  at  virtually  all  wave¬ 
lengths.  In  addition,  the  slowing  down  of  the  primary  electrons  produces  bremsstrahlung 
X-ray  emission.  The  calculation  of  ionization  produced  by  the  precipitating  electrons  is 
done  by  determining  their  spectral  distribution  from  the  measured  emissions.  Because 
we  are  here  restricting  our  discussion  to  altitudes  between  80  and  200  km,  the  electron 
energy  range  of  interest  is  approximately  500  eV  to  100  keV.  Once  the  electron  spectral 
distribution  is  known,  the  resulting  ionization  can  be  accurately  modeled  (see,  for  example. 
Vondrak  and  Robinson,  1985).  Thus,  remote  sensing  depends  on  the  accuracy  with  which 
the  energy  distribution  of  the  precipitating  electrons  can  be  determined  from  the  measured 
emissions. 


The  flux  of  X-rays  Fx  of  energy  Ex  is  given  by 


F.(Eg) 


${EX,  Et)Fe(Et)dEe 


(1) 


where  Ft  is  the  flux  of  electrons  of  energy  Ee.  The  function  $(Ex,Ee)  is  a  source  func¬ 
tion  relating  the  two  fluxes.  In  this  expression,  absorption  of  X-rays  in  the  atmosphere  is 
neglected;  this  is  justified  for  most  auroral  electron  spectra  where  the  average  energy  of 
the  electrons  is  well  below  about  40  keV.  Determination  of  the  electron  fluxes  from  the 
measured  X-ray  fluxes  involves  deconvolution  of  this  integral  equation.  There  have  been 
several  approaches  to  this  problem.  For  the  experiments  described  below  we  use  the  inver¬ 
sion  technique  described  by  Brown  (1972).  In  this  method,  the  differential  cross-section  for 
X-ray  bremsstrahlung  production  is  given  by  the  Bethe-Heitler  formula.  The  expression 
for  calculating  the  flux  of  X-rays  produced  by  a  given  distribution  of  electrons  is  in  a  form 
that  can  be  inverted  analytically.  Thus,  if  the  X-ray  fluxes  are  known,  the  electron  spec¬ 
trum  can  be  derived  from  the  measured  X-ray  spectrum  by  a  direct  calculation.  Because 
the  method  involves  taking  derivatives  and  second  derivatives  of  the  X-ray  distribution,  it 
is  sensitive  to  uncertainties  in  the  measured  X-ray  intensities.  In  the  data  used  here,  the 
measured  X-ray  fluxes  were  integrated  over  times  long  enough  to  ensure  good  counting 
statistics.  All  methods  for  deconvolving  X-ray  spectra  result  in  electron  spectra  with  large 
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uncertainties  at  the  lower  energies.  This  is  because  the  low  energy  X-rays  are  produced 
primarily  by  higher  energy  electrons.  Fortunately,  the  resulting  electron  density  profiles  in 
the  ionosphere  below  200  km  altitude  are  relatively  insensitive  to  the  low  energy  electron 
fluxes,  so  that  the  densities  can  be  determined  with  good  precision. 
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DESCRIPTION  OF  EXPERIMENTS 


The  ground-based  ionization  measurements  were  made  with  the  Chatanika  incoherent- 
scatter  radar  (Leadabrand  et  al.,  1972;  Baron,  1977;  Vondrak,  1983).  For  these  experiments 
the  radar  was  generally  operated  in  a  variety  of  modes  combining  elevation  scans  and  fixed- 
position  measurements.  The  ionization  was  measured  with  a  range  resolution  of  4.5  km 
and  with  integration  times  of  about  15  sec.  Because  the  satellite  and  radar  measurements 
are  generally  not  spatially  coincident  nor  simultaneous,  ground-based  all-sky  photographs 
were  used  when  possible  to  relate  the  measurements  made  during  the  nighttime  passes. 
Visible  auroral  imagery  from  the  DMSP-F2  satellite  was  not  available  for  the  cases  used 
in  this  study. 


X-ray  measurements  were  made  with  an  auroral  X-ray  sensor  on  the  DMSP-F2  satellite 
(Mizera  et  al.,  1978).  DMSP-F2  is  in  a  nearly  sun-synchronous  circular  polar  orbit  at  an 
altitude  of  830  km.  The  sensor  field  of  view  is  2°  by  14°,  providing  an  in-track  spatial 
resolution  of  about  50  km  and  a  cross-track  resolution  of  about  350  km.  This  corresponds 
to  about  0.5°  of  latitude  by  about  0.5  hour  of  local  time.  The  downlooking  X-ray  sensor 
provides  integral  X-ray  measurements  in  fifteen  energy  channels  from  1.4  to  20  keV.  These 
X-ray  spectra  were  converted  to  equivalent  incident  electron  spectra  using  the  inversion 
technique  described  by  Brown  (1972).  From  these  electron  spectra  we  calculated  the 
ionization  production  rate  q  as  a  function  of  altitude  using  the  method  of  Rees  (1963).  The 
electron  density  n  was  calculated  assuming  the  model  effective  recombination  coefficient 
a  given  by  Vickrey  et  al.  (1982).  A  transport-free  equilibrium  solution  was  used  for  the 
electron  continuity  equation,  such  that  the  electron  density  as  a  function  of  altitude  z  is 
given  by 


n(z)  = 


?(*) 


(2) 


M*)J 

For  the  sunlit  cases  the  production  rate  due  to  electron  precipitation  was  added  to  the  solar 
photoionization  rate  specified  in  the  empirical  model  of  Robinson  and  Vondrak  (1984). 


Precipitating  electron  fluxes  in  eight  differential  energy  channels  from  1-20  keV  were 
measured  once  per  second  by  the  J/3  instrument  on  DMSP.  Whenever  available,  the  pre¬ 
cipitating  electron  data  were  used  as  an  independent  test  of  the  validity  of  the  X-ray 
spectral  inversion  technique  (see  also,  Mizera  and  Gorney,  1982). 

The  entire  set  of  simultaneous  Chatanika/DMSP-F2  observations  was  searched  for 
cases  that  satisfy  the  criteria  of  spatial  and  temporal  coincidence,  with  aurora  present 


8 


within  the  Chatanika  field  of  view.  Another  constraint  was  that  the  auroral  luminosity 
and  ionization  had  to  be  fairly  homogenous  over  the  X-ray  sensor  field  of  view.  Only  four 
cases  were  identified  that  met  these  criteria.  In  all  cases  there  was  good  spatial  coincidence 
between  the  regions  of  enhanced  X-ray  luminosity  and  the  regions  of  enhanced  ionization. 
During  two  cases  the  conditions  of  simultaneity  and  homogeneity  were  sufficiently  good 
to  allow  a  detailed  quantitative  comparison  between  the  measured  and  inferred  altitude 
profiles  of  ionization.  In  the  next  section  we  present  in  detail  the  observations  made  on 
these  two  days  that  are  representative  of  auroral  conditions  in  the  evening  sector  and  in 
the  morning  sunlit  sector. 

OBSERVATIONS  OF  AN  EVENING  AURORA 

On  21  February  1979  the  DMSP-F2  satellite  passed  slightly  to  the  east  of  Chatanika 
at  0727  UT  (2127  local  time).  The  location  of  the  satellite  ground  track  is  shown  in  Figure 
1.  During  that  night  the  Chatanika  radar  was  operated  in  alternating  elevation  scans  and 
fixed-position  measurements.  The  E-region  limits  of  the  elevation  scan  closest  in  time  to 
the  satellite  pass  are  also  shown  in  Figure  1. 

During  the  DMSP  pass  there  were  bright  auroral  bands  south  and  overhead  of 
Chatanika.  Following  the  satellite  pass,  the  bright  auroral  bands  moved  southward.  The 
Chatanika  elevation  scan  began  in  the  south  within  one  of  the  bright  auroral  bands.  All-skv 
camera  photographs  taken  at  that  time  showed  that  the  auroral  bands  extended  over  the 
southern  sky  from  the  northwest  to  the  southeast  as  shown  by  the  shaded  region  in  Figure 
1.  On  the  figure  are  shown  the  times  corresponding  to  the  beginning  and  end  of  the  radar 
elevation  scan.  As  the  radar  scanned  northward  the  auroral  conditions  changed  dramat¬ 
ically  with  the  sudden  appearance  of  rayed  arcs  and  reduced  auroral  luminosity.  These 
temporal  variations  were  evident  both  in  the  all-sky  photographs  and  in  the  Chatanika 
measurements  made  in  the  fixed-position  mode  and  during  the  elevation  scan. 

The  radar  ionization  measurements  during  the  elevation  scan  closest  in  time  to  the 
satellite  pass  are  shown  in  Figure  2.  The  contours  show  electron  density  as  a  function  of 
altitude  and  invariant  latitude  in  the  magnetic  meridian  plane.  Conspicuous  features  are 
the  enhanced  ionization  between  62°  and  66°  latitude  that  is  associated  with  the  bright 
visible  auroral  band  and  the  high-altitude  ionization  in  the  north  between  66.6°  and  67. 2°. 
This  latter  ionization  is  associated  with  the  rayed  arcs  that  appeared  after  the  satellite  pass. 
The  peak  E-region  electron  density  in  the  bright  auroral  band  varied  between  3  x  105  cm~ 3 
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Figure  1.  DMSP  groundtrack  over  Alaska  on  21  February  1979.  The  location  of  the 
bright  auroral  band  at  0732  UT  is  shaded.  Also  shown  are  the  E-region  limits  of  the 
Chatanika  radar  elevation  scan  and  the  field  of  view  of  the  DMSP  X-ray  sensor. 
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and  6  x  lOscm  3  and  was  situated  at  altitudes  between  110  and  120  km.  Electron  densities 
in  the  rayed  arcs  maximized  at  altitudes  near  140  km. 

The  X-ray  sensor  on  the  DMSP  satellite  detected  enhanced  X-ray  fluxes  at  the  lati¬ 
tudes  near  Chatanika.  Because  of  the  temporal  variations  in  the  aurora  after  the  satellite 
pass,  we  compare  the  radar  and  satellite  measurements  made  in  the  bright  auroral  band 
to  the  south  of  Chatanika  'which  remained  relatively  stable.  The  energy  distribution  of 
bremsstrahlung  X-rays  in  the  auroral  band  south  of  Chatanika  is  shown  in  the  first  panel 
of  Figure  3.  The  measurements  are  shown  by  the  individual  points  with  error  bars  derived 
from  counting  statistics.  The  instantaneous  field  of  view  of  the  detector  is  shown  by  the 
rectangle  in  Figure  1.  The  X-ray  data  were  integrated  for  only  three  seconds  so  that  the 
total  field  of  view  was  about  40  per  cent  larger  in  the  direction  of  satellite  motion  than 
indicated  on  the  figure.  Because  of  the  large  field  of  view  of  the  X-ray  sensor,  the  measured 
X-ray  energy  distribution  represents  a  spatial  average  of  the  X-ray  fluxes  produced  in  the 
band. 

The  precipitating  electron  energy  distribution  derived  from  the  X-ray  distribution  in 
the  auroral  band  is  shown  in  the  middle  panel  of  Figure  3.  For  comparison  we  show  in 
figure  3  the  direct  measurement  of  the  precipitating  electron  flux  by  the  J/3  electron 
detector  on  DMSP.  The  J /3  electron  measurements  were  made  on  the  same  field  lines  as 
the  X-ray  sensor  field  of  view  and  were  averaged  over  a  time  interval  that  corresponds 
to  the  X-ray  sensor  field  of  view  along  the  satellite  ground  track  divided  by  the  satellite 
velocity.  Good  agreement  is  evident;  both  the  measured  and  derived  spectra  indicate 
average  electron  energies  of  about  7  keV  with  a  total  energy  flux  of  10.8  ergs/cm2-s. 

The  third  panel  in  Figure  3  shows  the  ionization  profile  computed  from  the  electron 
fluxes  that  were  determined  from  the  X-ray  data.  For  comparison  we  show  the  average 
electron  density  measured  by  the  radar  over  the  latitude  range  that  overlaps  with  the  X-ray 
sensor  field  of  view.  The  horizontal  bars  show  one  standard  deviation  from  the  mean  value. 
It  is  apparent  that  the  electron  density  was  fairly  uniform  within  the  X-ray  sensor  field 
of  view.  The  agreement  between  the  measured  electron  densities  and  the  derived  electron 
densities  is  very  good  at  all  altitudes.  We  expect  the  densities  computed  from  the  X-rays 
to  be  slightly  less  than  the  measured  densities  because  the  auroral  bands  did  not  quite  fill 
the  detector  field  of  view,  as  shown  in  Figure  1. 
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electron  flux  measurements.  The  horizontal  bars  on  the  radar  measurements  indicate 
the  results  of  spatially  averaging  the  data  over  latitudes  appropriate  to  the  X-ray 
detector  field  of  view. 


OBSERVATIONS  OF  A  SUNLIT  MORNING  AURORA 

On  15  June  1978  the  DMSP  satellite  passed  to  the  east  of  Chatanika,  intersecting 
the  Chatanika  magnetic  meridian  at  1547  UT  (0547  local  time  on  16  June).  The  DMSP 
ground  track  is  shown  in  Figure  4.  Because  the  day  of  this  experiment  is  close  to  the 
summer  solstice,  sunrise  at  Chatanika  was  more  than  four  hours  prior  to  the  DMSP  pass. 
The  solar  zenith  angle  at  Chatanika  was  about  70°  during  the  satellite  pass. 

The  Chatanika  radar  was  operated  in  continuous  sequential  elevation  scans  for  several 
hours  before  and  after  the  satellite  crossing.  Figure  5  shows  the  latitudinal  distribution 
of  ionization  measured  during  the  elevation  scan  made  at  the  time  of  the  pass.  A  broad 
region  of  diffuse  aurora  is  visible  over  Chatanika  with  maximum  ionization  of  about  3  x 
10s cm-3  between  115  and  120  km  altitude.  Thus,  the  average  energy  of  the  electrons 
producing  the  ionization  was  similar  to  that  observed  during  the  evening  aurora  described 
above.  Because  no  visible  auroral  imagery  was  available  for  this  sunlit  case,  we  examined 
data  from  successive  elevation  scans  to  determine  the  stability  of  the  auroral  precipitation 
during  the  satellite  pass.  In  Figure  6  we  show  the  maximum  E  region  electron  density 
measured  by  the  radar  during  nine  consecutive  elevation  scans.  The  data  are  plotted  as 
a  function  of  invariant  latitude  and  magnetic  local  time.  The  DMSP  trajectory  is  shown 
by  the  solid  line  slanting  toward  the  upper  left  part  of  the  figure.  The  location  of  the 
Chatanika  measurements  made  during  the  meridian  scan  that  was  simultaneous  with  the 
satellite  ground  track  is  shown  by  the  solid  line  in  the  left  part  of  the  figure.  Along  the 
satellite  ground  track  we  show’  a  rectangle  indicating  the  X-ray  sensor  field  of  view  for 
the  integration  period  during  which  the  maximum  X-ray  flux  was  measured.  The  dashed 
portion  of  the  rectangle  represents  the  increase  in  the  size  of  the  field  of  view  owing  to  the 
satellite  motion  during  the  six  seconds  over  which  the  X-ray  data  were  integrated.  Note 
that  the  field  of  view  does  not  overlap  the  Chatanika  radar  field  of  view  during  the  scan 
that  was  simultaneous  with  the  satellite  pass.  However,  the  data  from  successive  scans 
shows  that  there  was  only  a  small  change  in  maximum  E  region  density  at  later  local  times. 

The  X-ray  spectral  distribution  measured  by  the  DMSP  instrument  during  the  satel¬ 
lites  passage  over  the  diffuse  aurora  is  shown  in  the  first  panel  of  Figure  7.  As  in  Figure 
3,  the  vertical  bars  show  the  uncertainty  in  the  measurements  and  the  solid  line  is  a  fit  to 
the  data.  Comparing  this  spectral  distribution  with  that  shown  in  Figure  3  shows  that  the 
spectral  shapes  are  very  similar,  with  the  morning  sector  fluxes  lower  by  about  a  factor  of 
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Figure  4.  DMSP  groundtrack  over  Alaska  on  15  June  1978.  The  Chatanika  field  of  view 
during  a  meridian  scan  that  was  simultaneous  with  the  satellite  pass  is  also  shown. 
The  shaded  region  indicates  the  latitudinal  extent  of  the  diffuse  aurora. 
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seven  from  those  measured  in  the  evening  sector  event.  The  calculated  electron  spectrum 
has  an  average  energy  of  about  6  keV  and  energy  flux  of  1.5  ergs/cm2-s.  The  middle  panel 
of  Figure  7  shows  the  production  rate  of  ionization  resulting  from  the  derived  electron 
fluxes.  The  peak  ionization  rate  occurs  near  115  km  altitude  and  falls  off  rapidly  at  al¬ 
titudes  above  the  peak.  However,  solar  illumination  must  also  be  taken  into  account  in 
estimating  the  observed  ionization.  The  rate  of  ionization  produced  by  solar  illumination 
for  a  solar  zenith  angle  of  70°  was  determined  using  the  results  of  Robinson  and  Vondrak 
(1984)  and  is  shown  in  the  middle  panel  of  Figure  7.  Although  the  ionization  produced 
near  the  E  region  peak  by  sunlight  is  small  compared  to  that  produced  by  precipitating 
electrons,  the  solar  contribution  is  dominant  at  higher  altitudes.  The  resulting  electron 
densities  are  shown  in  the  third  panel  of  Figure  7.  Because  the  X-ray  detector  field  of  view 
encompassed  almost  the  entire  diffuse  aurora  at  the  time  of  the  measurement,  we  show 
for  comparison  the  range  in  electron  densities  measured  by  the  radar  across  the  diffuse 
auroral  region.  The  derived  profile  falls  well  within  the  range  observed  by  the  radar.  We 
note,  in  particular,  the  importance  of  adding  in  the  contribution  from  solar  illumination 
in  obtaining  the  agreement  at  the  higher  altitudes. 
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SUMMARY  AND  DISCUSSION 


The  comparative  analysis  of  the  Chatanika/DMSP  data  sets  shows  that  satellite-borne 
X-ray  measurements  can  be  used  to  determine  quantitatively  the  spatial  distribution  of 
electron  energy  deposition  and  ionization  in  the  auroral  ionosphere.  For  all  the  cases 
that  were  analyzed,  the  X-ray  measurements  could  be  used  to  identify  the  regions  of 
electron  energy  deposition  measured  independently  with  the  Chatanika  radar  and  the  all¬ 
sky  camera.  For  two  of  the  cases,  auroral  conditions  were  sufficiently  uniform  md  stable  to 
allow  a  detailed  comparison  between  electron  densities  determined  from  derived  energetic 
electron  fluxes  and  electron  densities  measured  directly  by  incoherent  scatter  radar.  The 
data  were  obtained  during  observations  of  an  active  evening  sector  aurora  and  a  morning 
sector  diffuse  aurora.  The  average  energy  of  the  precipitating  electrons  in  both  cases  was 
about  6  keV  in  the  most  intense  regions.  The  energy  flux  in  the  morning  sector  event  was 
lower  by  about  a  factor  of  seven  from  that  in  the  evening  sector  event.  This  difference  in 
the  calculated  energy  flux  was  consistent  with  the  factor  of  two  to  three  difference  in  the 
observed  E-region  peak  electron  density.  Because  the  morning  sector  aurora  was  sunlit, 
the  contribution  to  the  total  ionization  produced  by  photoionization  had  to  be  added  to 
bring  the  derived  and  measured  profiles  into  agreement.  This  is  the  first  time  quantitative 
remote  sensing  of  ionization  in  a  sunlit  aurora  has  been  accomplished. 

In  the  two  cases  described  in  detail  above,  the  X-ray  measurements  represented  spatial 
averages  because  of  the  large  field  of  view  of  the  detector.  Thus,  the  computed  profiles  were 
compared  with  electron  densities  measured  by  the  radar  during  a  scan  over  a  latitudinal 
range  consistent  with  the  X-ray  detector  field  of  view.  The  agreement  between  the  derived 
and  measured  densities  depends  on  a  combination  of  statistical  accuracy  of  the  measure¬ 
ments,  temporal  variations  of  the  ionosphere,  inhomogeneities  over  the  satellite  sensor 
field  of  view,  and  uncertainties  in  the  calculation  of  ionization  from  the  bremsstrahlung 
measurements.  An  additional  uncertainty  in  the  X-ray  remote  sensing  technique  is  the 
presence  of  proton  precipitation.  Although  electrons  axe  the  dominant  source  of  ionization 
in  the  aurora  on  a  global  basis,  precipitating  protons  can  occasionally  produce  comparable 
or  even  greater  amounts  of  ionization  (Basu  et  a/.,  1987).  Because  precipitating  protons 
produce  no  bremsstrahlung  X-rays,  their  presence  and  the  associated  ionization  cannot  be 
detected  by  the  DMSP  X-ray  imager. 

Despite  these  uncertainties,  it  is  obvious  from  our  results  that  the  X-ray  remote  sens- 
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ing  technique  holds  great  promise  for  monitoring  auroral  electron  precipitation  on  a  global 
basis.  Because  the  X-ray  technique  also  yields  information  about  the  spectral  distribu¬ 
tion  of  the  auroral  electrons,  these  measurements  can  be  used  to  determine  the  height 
profile  of  electron  density  at  altitudes  below  200  km.  Further  observations  of  the  type 
described  above  are  important  in  establishing  the  range  in  auroral  conditions  under  which 
the  technique  is  applicable. 
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